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1. Introduction and summary 

One of the deep mysteries of perturbative string theory is the nature of the Hagedorn 
transition. In all known string theories with spacetime dimensions greater than two, there 
exists a critical Hagedorn temperature, Th, above which the thermal partition function 
for free strings diverges. The divergence can be attributed to the exponential growth of 
the number of perturbative string states at high energies. It is generally believed that 
rather than a limiting temperature for string theory, the Hagedorn temperature signals a 
transition to a new phase [^^,0] , analogous to the deconfinement transition in a nonabelian 
gauge theory. From the worldsheet point of view, the Hagedorn temperature is associated 
with the appearance of new relevant operators in the worldsheet conformal field theory. 
String theory at a finite temperature can be described by strings propagating in a Euclidean 
target space with time direction periodically identified with a period given by the 
inverse temperature. The winding modes in the Euclidean time direction are irrelevant 
operators when the temperature T is small. The lowest winding modes become marginal 
at Th and become relevant above Th- The emergence of these new relevant operators 
signals that the worldsheet theory becomes unstable and most likely will flow to a new 
infrared flxed point. 

By studying the effective fleld theory near the Hagedorn temperature, Atick and 
Witten 1^ argued that the Hagedorn transition is a flrst order transition that happens 
below the Hagedorn temperature seen in perturbation theory. In particular, the genus 
zero contribution to the free energy is nonzero above the transition, i.e. the free energy is 
of order g~'^, where Qs is the string coupling constant. This implies that one can no longer 
ignore the back reaction of the thermal energy density on the background geometry, since 
no matter how small Qs is, the back reaction is always at least of order one. 

One expects that a precise understanding of the Hagedorn transition and physics of the 
high temperature phase should give us important insights into the fundamental structure 
of string theory. Given our very limited understanding of these important questions, it is 
thus of great interest to flnd a simpler setting that one could study similar questions in 
full detail. 

Recently, it was found that free A/" = 4 Yang-Mills theory with gauge group SU{N) 
on a three sphere also exhibits a Hagedorn type transition in the limit N oo 



^ See [0J8|,^,1O,11| for other recent discussions of phase transitions in weakly coupled Yang- 
Mills theory. 
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In the large N limit the number of gauge invariant operators grows exponentially with 
conformal dimension. This leads to a Hagedorn temperature Th, given by an order one 
constant times the curvature of S^, above which the thermal partition function diverges 
in the N = oo limit. At large but finite N, one finds a weakly first order phase transition 
in the perturbation theory, with the free energy of order 0(1) in the low temperature 
phase, while of 0{N'^) in the high temperature phase. 

From the AdS/CFT correspondence |[T^ , |T3| , [T^ , the Af = 4 SYM theory with gauge 



group SU{N) and coupling constant qym on describes type IIB superstring theory in 
AdS^ X S5. The string scale Ig and ten-dimensional Newton's constant G are given from 
Yang-Mills parameters asS 

where R is the curvature radius of the AdSs. The t'Hooft expansion in Yang-Mills 
theory corresponds to the Qs expansion (quantum gravitational corrections) in AdS, and a 
small t'Hooft coupling 9ym-^ implies a strongly coupled worldsheet. Free SU{N) Yang- 
Mills theory with large but finite appears to describe a string theory in AdSs x 5*5 in 
the limit 

gs^O, Is -^00, finite cx^ (1.1) 

i.e. we are working in the regime 

= 00 > i? > . (1.2) 

While very little is known about this theoryS, if we assume the validity of the AdS/CFT 
correspondence in this extreme case, we conclude that it should have the following prop- 
erties: 

• The theory has a perturbative stringy expansion in terms of genus expansion with an 
effective expansion parameter G/R^. 

• The theory has a Hagedorn divergence in perturbation theory due to exponential 
growth of perturbative stringy states. 

• The free energy is 0(1) below T^, while of order G~^ above Th- 



^ We omit order one numerical constants. 

It has been conjectured that this theory might be related to a theory of higher spins 

Sin,!!- 



2 



Even though this AdS string theoryi seems rather different from our famihar flat space 
string theories, the features listed above about the Hagedorn transition are qualitatively 
similar to those in flat spacetime. Here we have the advantage that the theory has a holo- 
graphic description in terms of a free Yang-Mills theory and thus its Hagedorn transition 
can be studied in detail non-perturbatively. 

The Hagedorn transition for free Yang-Mills theory is a large phase transition . 
The number of states at a given energy E can be enumerated by counting gauge invariant 
operators of conformal dimension E. Since the theory is free, this can be done exactly. At 
energies^ 1 <^ i? ^ A^, one can treat all operators of given dimensions as independent. 

E 

One finds that the number of states grows as e^H with Th given by an order one constant. 
At energies E ~ A^, the counting becomes more complicated, since finite A effects like 
trace relations between single and multi-trace operators are important. Treating all the 
gauge invariant operators as independent dramatically overcounts the number of states, 
resulting in the Hagedorn divergence. At energies E ^ A^, instead of counting gauge in- 
variant operators, it is simpler to treat the system as A^ species of gluons and quarks. Then 
standard results for the ideal gas implies that the number of states should be proportional 

1 3 

to e'^^ 2 E 4 ^ which is much slower than the Hagedorn growth. The Hagedorn transition can 
thus be understood as the crossover from a region in which the gauge invariant operator 
description is more appropriate to a region in which the gluon-quark description is more 
appropriate. On a compact space like S^, a sharp transition can only happen in the strict 
A = oo limit. At large but finite A, the transition arises as an artifact of the large A 
expansion0. The interpolation between the low and high temperature regimes should be 
completely smooth at finite A and becomes narrower and sharper as A is increased. 

Translating the above physical picture to the corresponding string theory in AdS, we 
conclude that 

• The Hagedorn divergence in the bulk string theory should be an artifact of the string 
perturbation theory. The transition should be smooth non-perturbatively. 

• The appearance of the Hagedorn divergence in perturbation theory can be attributed 
to a "stringy exclusion principle" , i.e. not all perturbative states are independent of 
one another at high energies. 



See also [17| for some recent discussions of Hagedorn transition in large AdS. 
^ We will take the radius of to be 1. 

^ More precisely, the partition function is an analytic function of temperature T for ah A^. The 
non-analyticity in T arises as a result of expansion in 1/N. 
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In this paper we study non-perturbative aspects of this Hagedorn transition using Yang- 
Mills theory. We would like to understand how non-perturbative stringy effects in AdS 
smooth out the sharp Hagedorn transition in perturbative string theory. In other words, 
we are interested in knowing how the Hagedorn divergence in free string theory is resolved 
at finite string coupling i. 

We now summarize the main results of the paper. We will argue that the large- 
Hagedorn transition for free A/" = 4 Super- Yang-Mills theory belongs to the same univer- 
sality class as that of the double-trace unitary matrix model 

Z{a) = J dUe^p [a{T)TrUTrU^ (1.3) 

where U is a unitary matrix and a(T) is a function of temperature. ( p..3|) has a Hagedorn 
type transition in the large limit at a(T) = 1. Thus the Hagedorn transition in the AdS 



string theory can be studied by examining the critical behavior of (|1.3| ) near a = 1. We 
study the matrix model ( p..3| ) to all orders in the 1/A^ expansion and find that there exists 
a double scaling limit 
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a- 1^0, A^^oo, (a - 1)A^3 = finite (1.4) 

in which the transition is smoothed out. Expressed in the language of the corresponding 
AdS string theory the limit is 

T-Tn^O, |^-.0, (T-TH)(^y= finite, ll^ = ^ (1.5) 

where l^p is the five-dimensional Planck length in AdSs. The physical meaning of ( |1.4| )- 
( p..5|) can be understood as follows. At infinite A^, there is a sharp phase transition at 
T = Th- At large but finite A^ (or R/l^p), the transition is smoothed out into a finite region 
around Th whose width scales inversely with A^ (or R/l^p) with some critical exponent. 
The above double scaling limit ensures that we stay within this region and decouples 
nonessential physics. 

In the double scaling limit (|1.4|) -( pr5| ), one finds a simple physical picture for the 
Hagedorn transition of the bulk string theory. This picture is reminiscent of the Hawking- 



Page transition |1^ ,19| in the strong coupling regime. We find that the Hagedorn divergence 



^ Naively gvM = seems to imply that string coupling is also zero. The effective genus 
expansion parameter gll^s/^ iii (1-1) for the bulk string theory is nonzero at finite N. This is 



the sense that we talk about the theory at finite string coupling throughout the paper. 
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is resolved by summing over two different classical string'?/ geometries. Since we are working 
in the Ig ^ oo limit (|1 . 2| ) for the bulk string theory, the concept of geometry in terms of 
classical gravity is not valid here. Nevertheless, in the large limit, classical stringy 
geometry, defined by exact worldsheet conformal field theory, is still a valid concept. 

We find that near the Hagedorn transition (both above and below), the full partition 
function of the theory can be written in a form 

Z = Zt + c-^^Zl . (1.6) 

Zt can be interpreted as the partition function for strings in an AdS background with 
imaginary time-direction periodically identified (we will call it thermal AdS), which con- 
tains a Hagedorn divergence. Z^ is the partition function of a noncritical string theory 
which emerges only as the Hagedorn temperature is approached. We find that it has 
the same scaling as the unconventional branch of the Af = 1 super-Liouville theory or 
a branched polymer. Sl can be interpreted as the difference in the classical string field 
action between the above two string backgrounds. When T < Th, e~^^ is exponentially 
small in large N and the thermal AdS dominates. At T ^ Th, Sl is of order one and the 
two string backgrounds are equally importantly. In particular, the Hagedorn divergence in 
Zt is cancelled by a similar divergence in Z^. In this regime, the theory does not have a 
(stringy) geometric interpretation. When T > Th^ e~^^ becomes exponentially large and 
the noncritical string background completely dominates. Zl is tachyonic below Th while 
the thermal AdS is tachyonic above Th- 

We also find some interesting non-renormalization properties for Zt and Z^. Zt 
contains only a genus one contribution with all higher loop contributions vanishing. Below 
Th, Zt also does not receive perturbative corrections beyond one- loop. Above Th, Zt 
does receive perturbative contributions to all loops. 

The plan of the paper is as follows. In section 2 we discuss the large Hagedorn 
transition of A/" = 4 SYM theory at infinite A^. In section 3 we discuss the Hagedorn 
transition at finite in Yang-Mills theory. We find that there exists a double scaling 
limit in which the partition function smoothly interpolates between the low and high 
temperature regimes. We discuss the interpretation of Yang-Mills calculation in terms 
of string theory in AdS. We conclude with discussions of implications of our findings in 
section 4. In Appendix A, we discuss a generalization of ( p..3|) which includes both single 
and double trace terms. This arises when one includes fundamental matter in the theory. 
We show that the double scaling limit of the model is the same as that of the single-trace 
unitary matrix model. 

^ More precisely, we mean T — Th ^ 0{N^^). 
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2. Free large Yang-Mills theory on 

In this section we briefly review the Hagedorn transition for a free SU{N) Yang- 
Mills theory (with adjoint matter) on in the inflnite limit While the results 
presented here are not new, our derivation of the phase transition is new. The basic idea 
is to integrate out all flelds in the theory except for the zero mode of the Polyakov loop. 
The partition function is then reduced to a unitary matrix integral. 

Expanding all flelds in Yang-Mills theory in terms of harmonics on S^, to the low- 
est order in qym, the theory reduces to a quantum mechanics problem of free harmonic 
oscillators 

C = lj2^r [{D.M^r - u^lM^] (2.1) 

a 

where Ma are N x N Hermitian matrices and the sum is over all fleld types and their 
Kaluza-Klein descendants on 5""^. iVa is the frequency for each mode. We will take the 
radius of 5""^ to be one. The covariant derivative in (^]l|) is given by 



DtMa = dtMa-i[A,Ma] 

A comes from the zero mode (i.e. the mode independent of coordinates on 5''^) of the 
time component of the gauge fleld and is not dynamical. Its equation of motion imposes 
the constraint that physical states must be SU{N) singlets, i.e. Gauss law on S^. The 
partition function of the theory at flnite temperature can then be obtained by integrating 



out all flelds in (|2.1| ) except for A. Projecting into the singlet sector of the standard results 



for harmonic oscillators one flnds 



Z= [duU {detadj (1 - eae-f'^^U)) (2.2) 

where U is an SU{N) matrix whose integration imposes the singlet condition^. It can 
also be interpreted as the Wilson line for A around the Euclidean time circle (Polyakov 
loop), detadj denotes the determinant in the adjoint representation and = 1 (—1) for 
bosonic (fermionic) Ma. Equation ( |2.2|) can be more conveniently written as 



Z = C{(3) j DUexp !^Trt/^Trf/-"j (2.3) 
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Note that we normalize the measure DU so that J DU = 1. 
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with 

Wr, = ZBin(3) + i-ir+hFin(3) . (2.4) 

zb{(3), zpif^) are the single particle partition function for the bosonic and fermionic sector 
respectively, e.g. 

Ml3)= Yl ^^(/3)= E • (2.5) 

a, bosons a,fermions 

The prefactor C{(3) in ( |2.3[ ) is given by 

oo 

logC{(3)=/3Eo-J2 — ■ (2-6) 
^ — ^ n 

n=l 

i?o is the zero-point energy of the theory on and the second term arises because the 
group is SU{N) rather than U{N). For free A/" = 4 SYM theory on 5"^, we have 

6a; + 12x^ — 2x^ 16x^ _n 

(1 - X)3 - "^=(1^' ^ = ^ • ^'-^^ 



The large limit of the matrix integral ( |2.3| ) can be evaluated using a saddle point 
approximation pil|PT]] . Depending on the values of Wn, it could have a complicated phase 
structure [^. For A/" = 4 SYM, one finds from ( p.7|) that parametrically 



^ « 1, n>l (2. 

Wi 



and only two phases arise [^J^] with the critical temperature Th given by the solution of 
the equation wi{Th) = 1. ( p. 81 ) implies that one can approximate ( |2.3D by the first term 
in the exponential 

[ DU e^-^'^^-'^' (2.9) 



In particular, since the critical behavior of ( |2.3| ) near = 1 is controlled by the first term, 
one expects that the two models ( p.3| ) and ( |2.9| ) should have the same critical behavior. 
Since we are interested in the transition region, for the rest of the paper we will focus on 



The matrix model ( p.9|) can be solved to all orders by introducing a Lagrange multiplier 
to eliminate the double trace term in the exponential^ 

Z{a)= I DU exp (aTtUTrU^) 

■ _ r 1 _ _ ,1 (2.10) 

DUdXdX exp 



1 



2na 



--XX + XTrU + XTrU^ 
a 



I thank M. Douglas and V. Kazakov for discussions on this point. A similar method was 
also used in the context of double trace deformation of hermitian matrix models by Klebanov and 



collaborators |23,24]. 
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AbsorbingEl the phase of A into U and letting 



lAI = -Ng 
II 2 ^ 



we find that 

Z(a) = , ^-^^ , 

2a Jo V 4a 

where F[g) is given by the unitary matrix integral 



Z{a) = — / gdge^v 7^ + N^F{g) ) (2.11) 



^N^F{g) ^ J DU exp (^^Ng{TtU + TrU^)^ . (2.12) 

Now Z{a) reduces to an integral of the partition function of a unitary matrix model 
weighted by a Gaussian factor over its coupling constant. 

The large expansion of matrix integral ( p.l2|) is well known [21,25, 26 1 and the 
leading order term is given by 



9' 



9<l 

F{g) = { - (2.13) 
ilog^-| g>l 

The discontinuity in the third derivative of F{g) at g = 1 is the Gross- Witten large N 
phase transition [ pT| , |27[] . The fi* < 1 phase corresponds to a saddle point of ( |2.12| ) at 



which the eigenvalues e*^% i = 1, ■ ■ ■ N ol U are distributed around the whole unit circle. 
In particular at g = 0, the distribution is uniform. In the g > 1 phase, the eigenvalue 
distribution develops a gap and does not cover the whole unit circle. When (7 oo all 
eigenvalues are localized at 6*^ = 0. 

The leading order result for Z{a) is now simply given by the saddle point of the g- 
integral ( |2.11j ). The Gaussian weight factor picks the eigenvalue distributions for different 



values of a from those of (|2.12|) . We introduce 

Q{a,g) = -^ + F{g) (2.14) 

and look for local maximums of Q. We find that: 
1. For a < 1, Q (a, (7) is a monotonically decreasing function of g with maximum at 

g = 0, Q(a,0)=0 (2.15) 
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This turns the SU (N) integral into an U (N) integral. 
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Expanding Q{a, g) around (7 = and performing the Gaussian integral we find that 



logZ(a) = -log(l-a) + --- (2.16) 
2. For a > 1, Q(a, 0) = is now a local minimum. The function has a maximum at 



^7o = :r^>l, w = Jl-- (2.17) 
1 — w \ a 

with 

TV^ / on \ 

log Z{a) = N'Qigo) = — + log(l -«;)+■•• (2.18) 

2 \1 — w J 

For AT = 4 SYM, a(T) is a monotonically increasing function of T 

a(T)=«;, = ^(^^, x = e-^ (2.19) 

with the critical temperature at 

w;i(/3h) = 1, /3h = 2 log(2 + V3) . (2.20) 

There is no contribution of 0{N'^) in the low temperature phase partition function ( |2.16| ) 
and the genus one term diverges as the critical temperature is approached, 

logZ^ -log(/3-/3H) + const, T^Th (2.21) 

It follows that the density of states around Th is (which can be found by a Laplace 
transform of ( |2.16| )) 

0(E) ^ const e^^-^ (1 + 0(1/^2)) (2.22) 



Note that ( p.21|) and ( p.22|) are typical of Hagedorn behavior. In fact they have exactly 



the same form as those in flat space string theory with all spatial directions compactified 
(see e.g. p8| , |29[] ). At order 0{N'^) there is a weakly first order transition at a = 1. The 



high temperature phase has nonzero genus zero contributionlli. 

In next section we will look at finite corrections to the above results. We would 



like to understand how finite corrections resolve the singularity in (2.21) and smooth 



Note that although ( ^.18| ) is continuous as a ^ 1, the second derivative of logZ(a) (i.e. 
specific heat) becomes singular. 
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out the transition. We conclude this section by making a few remarks on the physical 
interpretation of U in (|2.3|). 
1. U can be understood as the temporal Wilson line for A around the time circle r 
(Polyakov loop) 

r.n/3 

U"" = Pexp \ i I Adr] (2.23) 








The theory has a Zn symmetry, 

U-^e^U, keZ (2.24) 

coming from gauge transformations which are periodic up to an element of the center 
Zj\f. In terms of the eigenvalues e*^% i = 1, ■ ■ ■ N of U the Z^r symmetry corresponds 
to discrete translations 6i ^ 6i + ^j^. In fact ( |27^ ) has a larger "accidental ^7(1) 
symmetry" , since it was obtained from a free theory. The correlation functions of the 
theory, however, do not have to respect this larger symmetry, since the measure DU 
for SU{N) is only invariant. 
2. In the large limit, it is convenient to introduce density of eigenvalues 

1 ^ 

p{0) = -J2^i(^-^^)^ -^<0<ir (2.25) 



i=l 



Then 



and ( pT3| ) can be written as 

Z = C [ [Dp] e-^[^] (2.26) 



where V{p) has the form 

CXD ^ 

V[p] = N'J2-^\Pr^\" (2-27) 



n 

n=l 



Pn can be interpreted in the dual AdS string theory as the n-th winding modes around 
the time circle. From (|2.27[ ) and ( ^.20| ) , pi become massless at the Hagedorn tempera- 



ture and tachyonic above Th- The Hagedorn divergence (|2.21|) is precisely due to that 
Pi becoming massless [|[|[ . In the low temperature phase, the eigenvalue distribution 
is given by that of ( |2.12| ) at (7 = 0, which is 

= I ■ 
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In the high temperature phase, p{9) develops a gap. 
3. The Zn symmetry (|2.24|) becomes an ^7(1) in the large limit, corresponding to 



+ a (2.28) 
Under the transformation ( |2.28| ), pn transform as 



It is tempting to interpret the ^-circle as the dual time circle in AdS string theory. 
However, here a' = oo, it is not clear how to define T-duality. The U{1) symmetry 
of the originial time circle r ^ r + e is not manifest in the matrix model. Due to the 
symmetry ( |2.28|) , given a saddle point solution p{9) of ( p.26|) , p{9 + a) is a solution for 



any a. Thus there are a continuous family of saddle points in the high temperature 
phase. 



3. Hagedorn transition at finite 

In this section we would like to understand how finite N corrections in Yang-Mills 
theory resolve the singular behavior in the large expansion and smooth out the tran- 
sition. From the AdS/CFT correspondence, this tells us how non-perturbative stringy 
effects smooth out the singular Hagedorn behavior in perturbation theory for the bulk 
string theory. 

For this purpose, we need to work out the sub-leading corrections to the matrix integral 
( p.lOP near the transition region (i.e. a ~ 1). This can be obtained by expanding around 
the saddle points of 



A^2 /-oo / jst2 2 

Z(a) = -l ,„,exp(-^+iV^FM 



gage 



(3.1) 







2a 

using the large expansion of the unitary matrix model 



.N^Fig) ^ j DUexp(^^Ng{TrU + TrU^)^ . (3.2) 



We will first review the results for ^ ^ 
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3.1. Unitary matrix model 

The asymptotic expansion of ( pi .21 ) in large is well known pB| , [2E[ ] . Depending on the 
value of g, the expansion can be divided into three different regions 

/ Ar2^2 -I o^^-C/_^ no 1 — i(l) 



N^F{g) 



^ g I J_p-2iv/(s) J_ /?( 

4 ^ 2-K Z^n=l TV"-' " 



(2) 



(3.3) 



where 



(3.4) 



Explicit expressions for various series F^^ are not important for our discussion 

below and some important features of ( p.3|) are: 

1. There are no perturbative corrections to F[g) for g <1 beyond the leading term. All 
corrections are non-perturbative in A^. 

2. Sub-leading terms in the expansions for 5* < 1 and (7 > 1 become singular as (7 — ;> 1, 

1 1 



i?(3) 



(3.5) 



{1-g)^' " {g-lf^' 

The asymptotic expansions in break down and are replaced by that in interme- 

2 

diate region |(7 — 1| ~ N~3. 
3. One can define a double scaling limit 



g = l-N~3t, N ^ 00 

(2) (2) 
which extracts the F^ ' term from the expansion. F^ satisfies the equation 

^Ff (t) = -/^(t) 

where /(t) in turn satisfies the Painleve II equation 

1 



-r\t) = f+tf{t) 



(3.6) 



(3.7) 



(3.^ 



(2') 

Higher order terms Fn can be shown to satisfy more complicated differential equa- 

f2) 

tions. When |t| ^ 1, Fq has the following asymptotic expansion 

*i _ i Ino-r-fl - 3 -i- 63 
6 8 ^1 T28F 1024*6 

1 _ 4V2 



.(2) 



27r 



e 3 



1 _|_ 35 _ 3745 _|_ 

sVitl" 384F i8432V2ti 



-t > 1 
t > 1 



(3.9) 
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(2) 

Note that Fq is a smooth function of t and interpolates smoothly between g > 1 
{t < 0) phase and g < 1 {t > 0) phase. 
4. It was argued in that Fq ' (t) describes the full partition function of the type OB 
theory in d = dimension, i.e. pure 2-d supergravity. The parameter t is proportional 
to the cosmological constant // in the super-Liouville interaction. 

3.2. Critical behavior around a ~ 1 and a double scaling limit 

We are interested in understanding the critical behavior of the matrix model near 
the transition point a ~ 1. We emphasize that the integral ( |3.1| ) is manifestly finite for all 
values of a. The Hagedorn divergence arises as a result of doing the large expansion. 
To see how the Hagedorn divergence is cancelled in the finite theory, it is convenient to 
split the integral in (J3.1D into 



/ dgge — a— U ^)^Zi + Z2 
Jo 

1 — a V / 



(3.10) 



with 

Z2 = ^f gdge^ 



(3.11) 



1 



where F{g) is defined for g < 1 



F = N^F{g) - ^-4^ (3.12) 



In ( |3.10| ) we separated the integration for g into g < 1 and g > 1 regions motivated from 
( ^.3| ). We further isolated from the g < 1 part of the integral the term jzr^, which gives 
rise to the Hagedorn divergence ( p.l6| ), and a non-perturbative term 



1 — a 4a 

which precisely cancels the divergence at a = 1 . This latter term is not seen in perturbation 
theory. One can also readily check that both Zi and Z2 are convergent integrals for all a. 
They can be evaluated by saddle point approximations. Notice that F defined in ( |3.12| ) is 



exponentially small except near g = 1 (see ( p.3| ) and (|3.4| )). When a < 1, one finds that Z2 
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is exponentially smaller than Zx and Zx is dominated by a saddle point which approaches 
1 as a — s> 1 _ . When a > 1 , Zi is exponentially smaller than Z2 and Z^ is dominated by the 
saddle point ( |2.17| ) which again approaches 1 as a It thus follows that the critical 

behaviors of the theory around a = 1 can be obtained by evaluating around the 

region (7 ~ 1. 

To find the scaling region, let us first look at the limit of a — > 1 from the high 
temperature side ( |2.18|) . As a — > we find that 



Although logZ'(a) is continuous as a — > its second derivative with respect to a (i.e. 
specific heat) becomes singular in the limit. The leading non-analytic term N'^{a — 1)^ in 
logZ'(a) suggests a scaling a — 1 ~ 0{N~^). This also follows from a — 1 ~ ((/o — 1)^ and 



go = l + Va-l + 0{a-l), 




(3.13) 



the scaling of the unitary model ( p.6|) . This motivates us to define a double scaling limit 



1 - a = N~3s, 



00 



(3.14) 



With a change of variable 



(3.15) 



the exponent in ( |3.1| ) can be written as 




(3.16) 



-Nh + -st + Fl^'\t) + 0{N-^) 



After some simple algebra, one finds that 




(3.17) 



where Bn can be obtained from Fn . The leading term Bo{s) is given by 




(3.18) 



14 



where 9{t) is the step function. From ( |3.18| ) and (|3.9|) Bo{s) is a completely smooth 
function of s including s = 0. Higher order terms Bn, n > are increasingly complicated 
and we have not looked at them in detail. Our procedure should guarantee that all of them 
are well defined and smooth functions of s. They are not relevant in the double scaling 
limit dSJiD . 

In the double scaling limit ( 3.14 ), keeping only the leading term in ( |3.17| ), we have 



(3.19) 



where we have introduced 



Bois) = Bois) - - . 

s 



(3.20) 



Note that as s ^ 0, So is singular 



Bo{s) + ■ 

s 



(3.21) 



This singularity precisely cancels with that of the first term in ( |3.19|) . For s < 0, B{s) can 
also be written as 

/•mo 



Bois) 



while for s > 



Bois) 



dte^ 



(2) 



1 



(3.22) 



(3.23) 



One can easily generalize the above discussions to unitary models with both single 
trace and double trace terms, e.g. 



Z(a, yu) = / dU exp 



aTrUTrU^ + -NipTvU + jn:tU^) 
2 



(3.24) 



where U is a, unitary matrix. (|3.24| ) is relevant when we include Nf matter fields in the 
fundamental representation of the gauge group with Nf /N to be finite0. One can show 
that the critical behaviors of (|3.24 ) are in fact the same as those of the single-trace unitary 
matrix models. This implies there is no Hagedorn type transition for ( J3.24D 0. Some 
details can be found in Appendix A. 
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The actual matrix model resulting by integrating out all fields except for the zero modes of the 



Wilson loop is more complicated. But again one can argue that (3.24) is enough for understanding 
the critical behavior of the theory. 
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3.3. AdS interpretation 

The partition function Z ( |3.19| ) has the following form 

Z = Zt + c-^^Zl (3.25) 

with 



(3.26) 



log Zt = — log s + const 
log Zl = log Bo{s) + const 
and 

SL = Nis = N'^il - a(T)) ^ N'^a{TH){T - Th) (3.27) 

Zt is due to the saddle at (7 = (uniform distribution of matrix eigenvalues) and Zl comes 
from Zl + Z2 at (7 ~ 1. 

We would like to interpret equation ( |3.25| ) in terms of contributions from two different 
classical stringy geometries in the corresponding AdS string theory. The free Yang-Mills 
theory we are working with corresponds to the Ig ^ 00 limit ( |1 . 1| ) for the bulk string 
theory. Thus the concept of geometry in terms of classical gravity is not valid here. 
Nevertheless, since we are working in the large limit, classical stringy geometry, as 
defined by conformally invariant sigma-model on the worldsheet, is still a valid concept. 
In particular, saddle points (local maxima) in the matrix model (|2.3|) should correspond to 



worldsheet conformal field theories. The saddle points of (p.3[) in turn coincide with those 
of ( p.l|) in the scaling region (|3.14|) we are interested in here. Thus we shall interpret ( |3.25| ) 



as summing over two classical string backgrounds. Sl can be interpreted as the difference 
in the classical string field action for two backgrounds and Zt, Zl are partition functions 
around each background. Clearly log Zt should describe string theory in the thermal 
AdS geometry, i.e. AdS with the time direction periodically identified. The precise string 
theory interpretation of Zl is not clear to us. One possible identification is that log Zl 
describes the unconventional branch of the super-Liouville theory with s identified as the 
cosmological constant. In the next subsection we will also mention other possibilities. 

The physical picture refiected from (|3.25| ) can be summarized as follows. We first 
divide our discussions into three regions and then comment on some general aspects. 
1. T < Th, s > but not too close to zero. In this region the Zt (thermal AdS) 
dominates and Zl is exponentially suppressed. Note that Zt does not contain per- 
turbative corrections beyond one loop. Strictly speaking, Zl does not correspond to 
a local maximum of the full matrix model ( p.l|) . Nevertheless, we find it natural to 

16 



interpret it as corresponding to a conformal field theory on the string theory side. 
This theory appears to be tachyonic due to the negative sign in equation ( |3.21|) . For 
s ^ 1 we can evaluate ( |3.23|) as an asymptotic series in 1/s. Surprisingly, we find 
that the asymptotic series terminate beyond the Gaussian integration 

~ s"^ 5 / 1\ 

log So = logs + C" + ZTT + I terms nonperturbative in - J (3.28) 

192 2 Y s J 

where C is a real constant. This implies that log does not receive perturbative 
corrections beyond one-loop. The imaginary term iix in ( |3.28| ) again indicates that the 
theory contains tachyonic modes. More precisely, it has the behavior of a complex 
tachyon in 0-Euclidean dimension. 

2. The transition region: s ~ 0(A^-3), i.e. T -Th 0{N-^). In this region, the two 
terms in ( p.25| ) are of comparable strength. Since two backgrounds contribute equally 
in this regime, the full string theory does not have a well defined (stringy) geometric 
interpretation. The 1/s Hagedorn divergence of Zt is cancelled by a similar divergent 
term in Zl (see ( p.21|) ). We remarked below ( |2.27|) that the Hagedorn divergence 
in Zt can be understood as pi becoming massless. What is the interpretation of 
the 1/s singularity for Zl7 There are two possible interpretations. One is that 
( p.21|) arises from the volume factor of the Liouville theory, i.e. the theory becomes 
effectively noncompact when s 00. A second interpretation is that Zl develops 
a 0-dimensional complex massless mode (let us call it ai) as s — 0, given that the 
theory is tachyonic for s > and not tachyonic for s < 0. If true, it means that for 
s > 0, o"! is tachyonic, while for s < 0, pi becomes tachyonic. 

3. T > Tff, s < but not too close to zero. The saddle at g ^ 1 becomes the local 
maximum of the full integral. Zl is exponentially large and dominates. = is no 
longer a local maximum and Zt is not physically meaningfullll. When — s ^ 1, ( |3.22| ) 
can be evaluated as asymptotic expansions in 1/s using (|3.9|) 



log^o = Li)i + C--log(-s)+ 3 + '^^^ ^^ +o(^^] . (3.29) 

^ 3 16 ^ 96(-s)i 256(-s)6 \{-s)- ' ^ ' 



2 



Treating (— s) * as the effective string coupling constant, log Zl contains perturbative 
corrections to all orders. 



15 
16 



Of course, the numerical constant multiplying the volume factor has to match (p. 2 
Since we have omitted terms in ( p.l7| ) much bigger than Zt- 
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4. Sl (|3.27|) , which we computed exphcitly using the dual Yang-Mills theory, can be 



interpreted to be the difference of the classical string field action of two backgrounds. 
Sl cannot be computed in a first quantized formalism with strings moving in a fixed 
background. In a second quantized formalism, like string field theory, it is natural to 
expect Sl to be computable. 
5. Our discussion above was restricted to the region T — Th ~ 0{N~i). Going beyond 
this region requires understanding the full matrix model (|2.3| ) , which is a rather non- 
trivil task. We have studied the truncated model ( |2.10| ) to all orders in in detail 
for general a. We find that for T < Th, there are no perturbative corrections to the 
partition function of the thermal AdS background beyond one-loop. We expect this 
to be true for the full matrix model (|2.3|), given the quadratic nature of the action. 
We also find that in ( |2.1C1| ) there are non-perturbative corrections of the form 

g-TVlogW ^3_3Q^ 

and 

e-^'^(^) (3.31) 

( p.31|) again suggests contribution of another geometry. However, away from the 
scaling region, ( |3.31|) is subdominant compared to ( |3.3C1|) , so it is not clear whether 



it has an unambiguous meaning. While it is tempting to interpret ( |3.3CI| ) as due to 
D-instantons, it is not clear to us how to understand the logA^ factor. 

3.4- The unconventional branch of super- Liouville theory 

It is a natural question to ask whether log So has an interpretation in terms of a non- 
critical string theory0. The scaling behavior of ( p.29| ) corresponds to a string susceptibility 
exponent given by 

1st = \ (3.32) 

Here we discuss two possible, perhaps related, interpretations. 

A well known phase of random matrix models which has the value of exponent ( |3.32| ) 



is a branched polymer phase (see e.g. [^). The branched polymer interpretation seems 



to fit with the picture that in the high temperature (deconfinement) phase, quarks and 
gluons are "liberated" , and the continuous Riemann surface description breaks down 
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This subsection grew out of discussions with I. Klebanov and J. Maldacena. 
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The scaling behavior ( 3.32 ) also appeared in the double trace deformations of Hermitian 
matrix models P^ , P^ , |55| . Based on Feynman diagrams generated by the double trace term, 
it was argued in that the matrix model describes a branched polymer phase, when the 
coefficient of the double trace term is sufficiently large. In our case, the story is less clear 
since we do not have a good geometric picture of how discretized worldsheets arise in a 
unitary matrix model. It is also not clear to us how to interpret the scaling behavior of 
( p.28|) for s > from this point of view. 

We will now offer an alternative interpretation of ( |3.29| ), again drawing inspiration 
from the analogous questions in the context of double trace deformations of Hermitian 



matrix models. In it was found that a new phase appears when the coefficient of 
the double trace term takes a particular value. Klebanov and Hashimoto p3| , pl| gave an 
interesting interpretation of this new phase. They observed that the scaling behaviors of 
the new critical points in the presence of double trace deformations can be explained by 
changing the branch of the Liouville dressing. Our discussion of the double scaling limit, 
and in particular, the expressions ( ^.221 ) and ( |3.23| ) are rather similar to those obtained 
there. Given that Fq in ( |3.22|) and ( |3.23|) is identified with the conventional branch of 
the Af = 1 super-Liouville theory , it seems natural to identify log Bq with the other 
branch of the super-Liouville theory. We will now check that the scaling in ( |3.29| ) is indeed 
consistent with this proposal with (— s) identified with the Liouville cosmological constant. 
In Appendix A we study a double trace deformation of the simplest single trace unitary 
model and show that its critical behavior is the same as the single trace case except for 
the case we discussed in previous subsections. 

The super-Liouville action can be written as (we take a' = 2) 



S 



1 

2^ 

— / 

2tv 



1 



+ -QR^ + i)diP + i^di; + ^^h^e^^"^ + 2i^Q}?tl;iPe^ 



(3.33) 



with the central charge given by 



and h satisfying 



CL = l^ 2Q2 



Q 



For pure supergravity we need cl = 10 which leads to 



(3.34) 



(3.35) 
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There are two solutions to (|3.34|) , 

6_ = b+ = V2 (3.36) 

v2 

where 6_ is the standard branch which satisfies the Seiberg bound. b+ branch does not 
correspond to a local operator. 

The dependence of ( |3.33| ) on /io can be obtained by taking 

$ ^ $ _ 1 In ^0 
and we find the free energy on a genus h surface is 



4^' (3.37) 



For the branch b- we have 



The scaling agrees with that of the first line of ( |3.9| ) if we identify fxo > with —t. For 
Ho < 0, Fh are identically zero since there are only non-perturbative corrections in the 
second line of (|3.9| ). For the other branch we have 

F^, 4^'-'''^ . (3.38) 

This agrees with ( |3.29|) if we identify — s oc /xq > 0. It seems natural to identify ( |3.28| ) 
with the same theory with fiQ < 0. However, the term in ( |3.28[ ) does not seem to agree 
with the scaling of ( |3.38| ) for h = 0. There are several possibilities for the disagreemen tJii. 
First notice that term is analytic in a — 1, so this term might not be universal and does 
not have a Liouville interpretation. The second is that the scaling ( ^.38] ) does not hold for 
Ho < due to certain subtleties. 

It would be very interesting to have a more complete and precise string theory iden- 
tification of (^) and (g). 



It was also pointed out to us by I. Klebanov that ( p.28| ) has the same scaling as the c = —2 
matrix model. 
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4. Discussion and conclusions 

In this paper we investigated non-perturbative aspects of tlie Hagedorn transition for 
IIB string theory in an anti-de Sitter spacetime in the limit that the string length goes 
to infinity. We find that as one approaches the Hagedorn temperature perturbative string 
theory breaks down and a noncompact Liouville direction appears to open upS. In the 
double scaling limit T — Th ~ 0{N~i), N oo, the full partition function can be 
written as a superposition of those from the thermal AdS and the Liouville background, 
weighed by their relative classical action. The Hagedorn singularity of perturbative string 
theory around thermal AdS is cancelled by a similar divergence in the Liouville theory 0. 
Our discussions here apply to a variety of Yang-Mills theories or their string theory dual 
which have the same critical behavior as the double trace unitary matrix model (|2.10| ) . 

The summing-over-geometry behavior we see here is rather similar to that of the 



Hawking-Page transition [p!8| , p!9[ in the strong coupling. There are some differences. One 
difference is that we are summing over conformal field theories instead of classical gravity 
backgrounds. Another important difference is that in our case at a given temperature 
only one background is stable. In the strong coupling limit, the Hawking-Page transition 
appears at a temperature Tc ~ much lower than the Hagedorn temperature at 

Th ~ O(A^). When extrapolated to the strong coupling, the phase transition here may 
become the Hawking-Page transition or describe the transition between a meta-stable 
thermal AdS to an AdS black hole at Th- 

It is somewhat suprising and interesting that near the transition point, the bulk theory 
can be described by a Liouville theory. This is reminiscent of the appearance of long throats 



in other singular CFTs like the conifold [p7| , |38| , p9 



We also found some interesting non-renormalization properties. At temperatures be- 
low the Hagedorn temperature, the string partition function around the thermal AdS 
background appears to receive only one-loop contribution. Equation (|3.28| ) also does not 
receive corrections beyond one loop. It would be interesting to understand these properties 
better. 

Given the universal presence of the Hagedorn behavior in perturbative string theories, 
the lessons we learned here might serve as a useful guide for probing non-perturbative 



In this section we will assume the Liouville interpretation. 

This is also reminiscent of tlie discussion in |36] where unitarity of the AdS black hole 



background is restored by summing over other geometries. 
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stringy effects in other string theories including those in asymptoticaUy flat spacetime. For 
example, it is interesting to check whether scaling behavior exists in type II or Heterotic 
theory in flat spacetime as one approaches the Hagedorn temperature, i.e. we would like 
to know whether leading order Hagedorn divergences at higher genera have the form 

oo ^ 

XogZ ^ -log(/3 - M + (4.1) 

with 7 the critical exponent. We hope to come back to this question later. 
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Appendix A. Mixed unitary matrix model 

Here we consider the double scaling limit of the following unitary matrix model 



Z{a, n) = I dU exp 



which arises if one introduces fundamental matter [0,^ . Our discussion follows that of ||24 
Introducing the Lagrange multiplier we flnd that 



8na 

/-oo 



dUdXdX exp 



iV2 

4a 



(A -fi){X--p) + [XTrU + XTrU^) 



Ana 



gdg 1^ d9 exp (^-^ {g' - g{^e-^' + Jle^') + + N'F{g)^ 



N 
2^ 



2 /"OO 



where 



gdge 



N'Q = log/o (^^) - ^(IH^ + g') + N'F{g) 



(A.l) 



and we have used 



2^'^(f) = j dU exp (^NgiTrU + Trt/^)^ 
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and 

de 



Note that we have introduced g = |A| and in the second hne of ( A.l ) we absorbed the 



phase of into the integration of 6. Below for notational simphcity we wiU denote |//| 
simply as ii. Thus ^ > 0. 
Note that 



We will evaluate (|A.1|) using the saddle point approximation. We will look at the 
leading order terms in Q. Using (|3.3| ) we find that 

1. When a < 1 and fj, < hq = 1 — a, Q has a maximum at 

/U ^/ N /"^ n'^ 1- a 

fio 4(1 - a) /Uq 4 

2. When a < 1 and /j, > hq or a > 1, Q has a maximum at 



2 



^o = a + | + W(a+ - a, Q(^o) = ^ (^o + - 1 - log^o) - ^ 



Thus it follows from the above that in the a — /i plane, below the line a + = 1 (note 
/i > 0), the system is in the phase in which the distribution of eigenvalues of U has no 
gap on the unit circle, while above the line a + fx = 1, the system develops a gap in the 
distribution of eigenvalues of U. 

The behavior of the system as the critical line is approached can be obtained as follows. 
It is clear that as far as /U 7^ 0, the saddle points of both phases approach (70 = 1 as the 
critical line is approached HI. More explicitly, from the above equations one can check that 
from both sides 

A 

^70 = 1 + :^ + 0(A2), A = ;U-^o<1 

1 — a 



Motivated from the scaling behavior of the unitary matrix model ( p.6| ), we let 
and 

N^Q = Ko-^'-^it- sf + (t) + 0{N-i) 
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This is different from /i = case, in wliich the saddle point is always at g = for a < 1. 
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with 

Ko = —{l-a) + ^ >- + —s\l - a) - - lo^ 

Plugging the above expressions into ( |A.1D we find that 



4 



We thus find that 



/OO 2 
-OO 



2 , 



log Z = + i^o^^ {s) +0{N-i\ 



with non-universal terms 

= ^(1 - a) + ^ ^ + " ~ - a) 

Thus we find that with any nonzero /i, the double trace deformation of the unitary matrix 
model has the same critical behavior as that of the unitary matrix model up to non- 
universal terms. The case for fi = needs separate treatment, which is discussed in the 
main text. 
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